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ROUGH  SURFACE  BACKSCATTER 


The  intent  of  this  contract  was  to  further  analyze 
some  model  computations  of  backscattering  strength  which  were 
done  in  summer  of  1987.  This  computational  model  * (equation 
(1)  below)  bracked  (with  a  3-5  dB  difference)  the  Chapman- 
Harris-Scott  empirical  model  for  backscattering  over  a  wide 
range  of  parameters,  (150-1200  Hz,  10-45  deg.  grazing,  and 
10-40  knots  winds) .  However,  the  computation  was  not 
properly  windowed,  and  when  a  suitable  window  was  introduced, 
the  computational  backscattering  strength  agreed  very  closely 
with  the  standard  perturbational  model  for  backscattering, 
(a  model  valid  only  for  small  values  of  the  Rayleigh 
parameter) . 


This  result  is  in  itself  rather  surprising,  as  the 
integral  involved  (derived  via  a  Kirchoff  approximation)  has 
the  form 


SS 


ka  qz 


2  ffd2„  e-'V’ftlP)  [e-02z°2<1-W<°)]-e-'5z2-’2l/(4*>2 


So 


(1) 


where  qz  =  2k  sin  (9),  k  =  acoustic  wave  number,  9  =  grazing 
angle  relative  to  mean  surface,  q±=  (-2kcos(e),  0),  P=  (px, 
Py)  ,  Sq  =  illuminated  area  of  surface,  A(p  )  =  windowing 

function  (e.g.,  raised  cosine,  etc.),  o2  =  surface  height 
variance,  and  W(p  )  =  normalized  surface  autocorrelation 

function.  The  computed  result  is  very  nearly  independent  of 
Sq  and  A ( p ) ,  for  suitably  tapered  windows, 

Clearly  for  o2qz2  <  <  1,  the  above  integral  reduces 
to  the  perturbation  result 

SSpert  ’ 
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(2) 


o2qz2  > 


where  F(K)  =  the  surface  power  spectrum.  But  for 
1,  it  is  by  no  means  obvious  why  SSka  s  SSpert'  (within  2-3 
dB)  .  Note  that  for  1200  Hz,  45  deg,  and  40  knots,  o2qz2  = 
50. 


As  is  well  known,  the  SSpgrt  scattering  model 
depends  only  weakly  on  frequency  and  wind  speed,  (using  the 
Toba  or  Donelan-Pierson  spectrum) .  This  is  in  contrast  to 
many  data  sets,  which  exhibit  a  fairly  strong  dependence  upon 
both  parameters.  Due  to  the  computational  results  cited 
above,  it  follows  that  the  SSka  model  has  this  same  weak 
dependence . 

It  should  be  noted  that  both  models  (1)  and  (2)  are 
incorrect  for  large  Rayleigh  parameter.  t^ie  first 
term  in  an  expansion  in  the  Rayleigh  parameter,  and  is  only 
valid  for  small  values,  where  the  higher  order  terms  are 
negligible . 

The  SSka  model,  and  in  fact  all  analytical  scatter¬ 
ing  models  that  use  integrals  of  ^Pbdry'/^n  (the  normal 
derivative  of  the  pressure  field  in  the  boundary) ,  are 
inaccurate,  because  these  models  must  assume  some  functional 
representation  for  dP^  /dn,  in  general  unknown.  Recall 
that  the  boundary  pressure  is  required  to  be  zero  on  the  free 
boundary.  The  Kirchoff  approximation  assumes  the  field  in 
the  boundary  is  that  due  to  the  source  and  an  image  source  in 
the  reflected  region.  But  this  field  is  only  zero  on  the 
flat  mean  surface,  not  on  the  actual  ocean  surface  itself. 
This  approximation  implies  that  ^^dry7^11  =  23P^nc/3n.  s^nce 
modifying  the  boundary  field  to  zero  may  well  require  large 
changes  in  the  derivatives,  there  is  no  way  of  assessing  the 
error  in  this  approximation. 
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The  Kirchoff  approximation  may  also  be  thought  of 
as  the  approximate  solution  of  an  integral  equation  satisfied 
by  ^bdry/^n,  of  the  form 


[I-T]  [3P 


bdry 


/bn] 


=  23P.  /bn 
me 


where  T  is  an  integral  operator  over  Sq. 

Now  [I-T]  1  =  2Tn,  0  <  n  <  ®  ,  if  ||T||  <  <  1,  over 
some  function  space.  If  the  Neumann  series  converges,  then 
the  Kirchoff  approximation  is  equivalent  to  the  approximation 
[I-T]  -  I.  But  for  large  Rayleigh  parameter,  i)  it  is  not 
known  if  the  series  converges;  and  ii)  no  error  bounds  are 
known  on  the  terms  Tn  ^Pbdry^n^  * 

With  respect  to  obtaining  a  second  order  theory, 
Winebrenner  in  his  thesis  uses  a  Rytov  method  of  partial 
summation  to  obtain  a  scattering  model  that  includes  second 
order  terms  in  the  phase  term.  But  Winebrenner  also  shows 
that  his  model  reduces  to  sspert  for  small  Rayleigh 
parameter,  and  to  SSka  for  large  values.  Since  the  computa¬ 
tional  model  shows  that  these  two  models  agree  for  the  ocean 
surface,  there  is  no  reason  to  believe  Winebrenner' s  model 
gives  anything  new  for  the  ocean  surface. 


Another  class  of  theoretical  scattering  models  are 
the  composite  surface  models,  SS  ,  where  the  scattering  is 

v->  O 

regarded  as  diffraction  from  a  small  surface  (high  wave- 
number)  ,  riding  on  a  large  surface  (low  wavenumber) . 


Since  one  rationale  for  introducing  SS  is  to  more 

C  o 

accurately  evaluate  the  integrals  in  SS^a,  the  recent  compu¬ 
tations  have  reduced  the  need  for  such  an  approach.  Other 

derivations  of  SS  are  ad  hoc,  and  do  not  begin  with  SS.  , 

cs  K.a 
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but  the  resulting  SS  model  agrees  with  those  which  begin 

C  5 

with  SSka. 

However,  if  one  does  introduce  a  large  and  a  small 
surface,  separated  by  some  dividing  wave  number  KD,  the 
result  is  that  in  a  number  of  the  standard  SS  (see  Brown, 
Bahar  and  Barrick,  McDaniels,  the  NRL  model,  etc.),  the 
dominant  effect  is  the  introduction  of  an  "effective"  grazing 
angle,  i.e.,  e  is  replaced  by  0  +  8  where  tan(6t^lt)  = 

the  rms  slope  of  the  large  surface. 

This  model  does  produce  more  wind  speed  dependence 
and  frequency  dependence,  depending  upon  how  one  chooses  the 
divider  KD,  but  not  enough  of  an  increase  to  adequately  model 
the  data. 


To  add  to  the  above,  a  new  composite  surface  model 

has  been  developed  by  the  San  Diego  group  (Dashen-details 

not  available) ,  which  chooses  KD  very  small,  so  that  the  tilt 

of  the  large  surface  is  negligible.  Of  course,  in  this  case, 

SS„„  -  SS_  .  once  more, 
cs  pert 

Surveying  all  the  existing  scattering  models,  one 
is  forced  to  the  conclusion  that  there  is  no  reason  to 
conclude  that  one  of  them  is  accurate  for  large  Rayleigh 
parameters . 

Due  to  the  above  difficulties  with  surface  scatter¬ 
ing  models,  and  also  the  problem  with  zero  Doppler  back- 
scattering,  there  has  been  a  recent  renewal  of  interest  by 
DARPA  in  bubble  scattering  models.  The  JASON  committee,  in 
the  summer  of  1987,  considered  and  rejected  a  bubble  model 
based  upon  scattering  from  bubbles  with  diameters  in  the 
micron  range. 
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In  the  spring  and  summer  of  1988,  consideration  was 
given  to  scattering  from  sausage  shaped  bubbles  (several 
centimeters  in  diameter)  which  detach  from  the  base  of 
capillary  waves  as  the  wind  increases,  (see  Toba,  Longuet- 
Higgins) .  No  conclusions  are  available  concerning  this 
model,  as  apparently  too  many  of  the  relevant  model 
parameters  are  unknown,  (e.g.,  density,  duration,  depth  of 
the  bubble  layer  below  the  surface,  etc.). 

More  recently,  attention  is  being  paid  to  a  third 
bubble  model,  consisting  of  large  foam  bubbles,  (meter 
diameter)  which  slide  down  the  face  of  large  waves  beneath 
the  surface.  As  of  this  writing,  the  results  of  this  bubble 
model  are  also  unknown,  again  most  probably  due  to  the  lack 
of  knowledge  of  the  governing  parameters.  Because  of  this 
uncertainty,  it  may  be  some  time  before  conclusions  are 
obtained. 

Even  if  ultimately  some  bubble  scattering  model 
successfully  explains  the  data  (most  particularly  zero 
Doppler  scattering) ,  an  accurate  surface  scatter  model  now 
would  certainly  clarify  the  need  for  additional  scattering 
mechanisms,  such  as  bubbles,  in  order  to  satisfactorily  model 
the  data. 


Accordingly,  during  the  last  stages  of  this  con¬ 
tract  a  new  computationally  intensive  approach  towards  an 
exact  scattering  model  was  begun.  Using  a  computer  simulated 
ocean  surface,  an  exact  (to  2-3  place  accuracy,  or  as 
required)  scattering  strength  for  that  surface  is  obtained. 
No  analytical  assumptions  regarding  are  made,  and 
an  error  criterion  is  developed  to  evaluate  the  accuracy  of 
the  computed  scattering  strength. 
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Performing  this  calculation  for  N  surfaces  (N  =  30, 
40?)  will  give  enough  data  to  estimate  the  ensemble  averaged 
scattering  strength  (the  usual  notion  of  scattering  strength, 
the  ensemble  average  of  the  2nd  moment  of  pressure)  .  This 
will  be  the  first  valid  estimate  of  scattering  strength  for 
large  Rayleigh  parameter. 

In  addition,  these  computations  will  also  provide 
estimates  of  the  variance  of  scattering  strength  (from  a 
theoretical  point  of  view,  this  requires  the  ensemble  average 
of  a  4th  moment  of  pressure) .  This  last  statistic  should  be 
a  valuable  addition  to  our  knowledge  of  the  surface  scatter¬ 
ing  problem. 

Major  portions  of  the  coding  for  this  approach  have 
been  completed,  and  preliminary  results  indicate  the 
theoretical  approach  is  correct,  and  the  computations  out¬ 
lined  above  are  feasible. 

Because  this  method  is  implemented  on  a  computer, 
it  is  flexible,  and,  for  example,  rather  than  using  a  finite 
Fourier  series  as  a  model  for  the  ocean  surface,  could  easily 
be  modified  to  include  in  the  surface  model  the  large  bubbles 
considered  in  the  third  bubble  model. 

-> It  also  seems  very  likely  that  this  approach  can  be 
modified  to  address  scattering  from  a  rigid  boundary,  i.e., 
the  rough  ocean  basement.  As  more  accurate  spectral  models 
of  the  ocean  basement  are  becoming  available,  computer 
simulated  basement  surfaces  can  more  accurately  model  the 
actual  surface.  And  the  above  computational  model,  extended 

>  *v>  .  •  -  . 

to  the  case  of  a  rigid  boundary,  could  provide  insight  into 
such  questions  as  the  anomalous  time  spreads  observed  in 
scattering  from  the  rough  basement. 

>  'Ti  ' 
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